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ABSTRACT 

I present a specific worked example of evolution through inflation to the initial con- 
ditions for an isocurvature CDM model for structure formation. The model invokes 
three scalar fields, one that drives power law infiation, one that survives to become the 
present-day CDM, and one that gives the CDM field a mass that slowly decreases during 
infiation and so "tilts" the primeval mass fiuctuation spectrum of the CDM. The func- 
tional forms for the potentials and the parameter values that lead to an observationally 
acceptable model for structure formation do not seem to be out of line with current ideas 
about the physics of the very early universe. I argue in an accompanying paper that the 
model offers an acceptable fit to main observational constraints. 

Subject headings: cosmology: theory — cosmology: large-scale structure of universe — 
galaxies: formation 
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1. Introduction 

This paper with its companion (Peebles 1998; here- 
after Paper II) is the latest in a series of studies of 
isocurvature models for structure formation. (Ear- 
lier papers may be traced back from Peebles 1997). 
There are two motivations for this work. First, it is 
important to know that there is an observationally 
viable and theoretically not unreasonable alternative 
to the commonly discussed adiabatic cold dark matter 
(ACDM) model for structure formation. As long as 
there are alternatives it demonstrates that we do not 
have an established standard model for the early uni- 
verse, that there is more to be discovered than tighter 
constraints on parameters and functional forms for a 
potential. Second, the isocurvature models seem to be 
better adapted to galaxy formation at high redshift, 
a condition I find attractive. This paper is meant to 
demonstrate that initial conditions for the isocurva- 
ture case, which I shall call ICDM, can be given an ex- 
plicit and not unreasonable basis in a physical model 
for inflation. In Paper II I argue that the model offers 
an acceptable fit to the observational constraints that 
can be applied without the use of numerical simula- 
tions. 

There is considerable activity in the study of the 
rich variety of ideas for a specific model for inflation 
motivated by current ideas in particle physics (as re- 
viewed by Randall 1997). I make use of functional 
forms for potentials that commonly appear in these 
discussions, but the motivation is to obtain a spe- 
cific model that could have come out of inflation to 
compare to the rich suite of observational evidence. 
The two cultures — based in particle physics and as- 
tronomical phenomenology — will meet if it becomes 
possible to establish that the narrow range of mod- 
els from fully acceptable physics overlaps the narrow 
range of models that are observationally acceptable. 

The dynamical actors in the structure formation 
model to be discussed here and in Paper II are the 
same as in the family of adiabatic CDM (ACDM) 
models — baryons, radiation (the CBR) with initially 
homogeneous entropy per baryon, cold dark matter, 
and three families of neutrinos — and in about the 
same amounts. In the adiabatic version, a scale- 
invariant departure from a homogeneous primeval 
mass distribution has power spectrum P (x k. In 
the ICDM model a scale-invariant spectrum of the 
distribution of the CDM would be P oc k~^, while 
the net mass density is homogeneous. Early discus- 



sions of inflation accepted the proposition that adi- 
abatic and isocurvature scenarios are equally well 
motivated (eg. Steinhardt & Turner 1983, Linde 
1985; Seckcl & Turner 1985), but when Efstathiou 
& Bond (1986) demonstrated that a scale-invariant 
isocurvature model violates the bound on the CBR 
anisotropy attention naturally turned to the scale- 
invariant adiabatic case. The COBE detection of the 
CBR anisotropy (Smoot et al. 1992) showed that 
if the universe were Einstein-de Sitter the ACDM 
spectrum with a reasonable bias would have to be 
"tilted" from scale-invariance, an arguably natural 
adjustment of the inflation picture (eg. Crittenden 
et al. 1993). Tilt is not needed if the mean mass 
density (excluding a term in the stress-energy tensor 
that acts like a cosmological constant) is well below 
the Einstein-de Sitter value, but the precedent has 
been set: consider tilting the ICDM spectrum to fit 
the observations. Models for inflation that tilt the 
spectrum in the wanted direction have been discussed 
by Kofman & Pogosyan (1988), Salopek, Bond, & 
Bardecn (1989), and Linde & Mukhanov (1997), and 
a model fitted to the observational constraints has 
been presented in Peebles (1997). Here and in Pa- 
per II I present a more detailed discussion along the 
lines of this last paper. 

For the purpose of displaying a specific worked ex- 
ample I adopt definite values of the parameters in the 
cosmology and the structure formation model and de- 
rive from them the parameters in the inflation model, 
the latter including the needed initial conditions early 
in inflation. The parameters will have to be reconsid- 
ered with each improvement of the observations, of 
course. The hope is that such adjustments driven by 
advances in the observations and physics may back us 
into the corner of model and parameter space that is 
a reasonable approximation to reality. 

The elements of the inflation scenario are presented 
in §2. Section 3 starts from adopted values of the pa- 
rameters in the models for cosmology and structure 
formation and prc^scuits derived values for the param- 
eters in the inflation model. Except where otherwise 
indicated units are chosen no fi = 1 = c, and I follow 
the notation in Peebles (1993), §17. 

2. The Inflation Model 

The model assumes three real single scalar flelds, 

one whose energy density drives homogeneous infla- 
tion, one that survives as the present-day CDM, and 
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one that serves to give the CDM field a rolhng mass 
that tilts the spectrum of quantum fluctuations frozen 
into the CDM distribution during inflation. A sim- 
ple way to tilt to a near power law power spectrum 
of frozen fluctuations starts with the condition that 
during inflation the universe is expanding as a power 
of time (Abbott & Wise, 1984; Lucchin & Matarrese 
1985), 



a oc t 



H ^ a/a^ [ety 



(1) 



The past event horizon is eliminated (or rather can 
be very large) if e < 1; a convenient value in the 
numerical example below is e = 0.1. 

The purpose of this discussion is to obtain a spe- 
cific model for the relevant properties of the cold dark 
matter. I do not address the issue of how the energy 
density in the field that is driving infiation is con- 
verted to the entropy that ends up mainly in the CBR 
today. 

2.1. Power Law Inflation 

Following Lucchin & Matarrese (1985), the power 
law expansion law in equation (^) assumes the energy 
density during inflation is dominated by the field x(t) 
with potential energy density that is an exponential 
function of the field, 



M2 



,-x/M 



Then the field equation is 

(fx 3 dx _ dV _ M 



-X/M 



(2) 



(3) 



for the power law expansion in equation (|l|), and the 
adopted solution is 

gX/(2Af) ^ (^/^^)(g_2e)-i/2. (4) 

The energy density in the field is 

UP 



(5) 



where the dot means derivative with respect to proper 
time t. Consistency of the Friedmann equation 



with equations (^ and (||) requires 

A-P _ 1 
mpi^ 167re' 



(6) 



(7) 



where mpi — G ^P is the Planck mass. Thus the 
mean mass density is 



3"^^i 
87re2t2 ■ 



(8) 



Given e, equation (Q) fixes the mass M in the po- 
tential Vy^ in equation (||) and sets the scale for x- 
Then the constant fixes the time scale for the evo- 
lution of X it passes through x = 0. 

2.2. The Cold Dark Matter 

If the Hubble parameter H is close to constant dur- 
ing infiation and the mass of the CDM field (p is much 
less than H then the spectrum of fiuctuations frozen 
into (f> is the scale- invariant form cx k^^ . To fit 
the observations the fluctuations in must be tilted 
to favor smaller scales (Peebles 1997). This is done 
by coupling to a field ^ that gives (j) a mass com- 
parable to H. The mass causes decay of fluctuations 
of (f> on scales larger the Hubble length, thus favoring 
the shorter wavelengths that pass through the Hubble 
length later. 

The Lagrangian density for the two fields is 



where the potential energy density is 



(9) 



(10) 



Here (5 and 7 are dimensionless constants; in the nu- 
merical example 7 ~ /3 e. As will be described, the 
quartic term in f/^ causes the effective mass of 4> to roll 
to lower values at a rate that can be adjusted to tilt 
the fluctuation spectrum of to an observationally 
acceptable power law, and the constant mass /x pro- 
duces a break in the power law to make (0^) converge. 
The constant mass m is assumed to be much smaller 
than iJ during inflation. After inflation, when -0 has 
dropped to zero, the CDM field has mass density 



p(x) = \m^4>{^'^ 



(11) 



where 0(x) is the field frozen (or squeezed) during 
infiation. This is the primeval mass distribution that 
seeds structure formation in the ICDM model. 

The potential in equation ( [To| ) is introduced ad hoc 
to solve the problem of tilting the CDM fluctuation 
spectrum, but it will be argued in §4 that it does have 
a pedigree in similar forms that appear in different 
contexts in many discussion of inflation. 
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It will be assumed that during inflation the field ^ 
is little affected by its coupling to (/), ■0 is a function 
of time alone, and m <^ tp. Then the wave equations 
from equations and ( ]lO| ) are 



(fib 3 dip ^ , o o , 



d'^(t) 3 d(t> 
W ^ Tt'dt 



With the changes of variables 

a — aiT^^'^, K — fc/(/iai), 



(12) 



(13) 



Equation (|17|) assumes space is periodic in some large 
coordinate volume Vu, so the creation and annihila- 
tion operators satisfy the discrete relation [at, aj^,] = 
5k,k'- With the normalization in equation ( p^ ) the 
field momentum vr — a^dip/dt satisfies the canonical 
commutation relation [(/)(x, i), 7r(x', t)] = i(5(x — x'). 

With the initial condition in equation ([l^), the 
mode amplitude 4>k{t) approaches the constant value 
(j)"^ at large r, and the field is frozen as a random 
Gaussian process with autocorrelation function 



^0(2^12) = ((^(xi)(^(xi)}, 
where the power spectrum is 



(19) 



where k is the comoving wave number of the Fourier 
mode (pkC^^'^, the field equations are 



3^* ^ . „ 
e-^ + -— + + = 0, 



dr2 



3 dcpk 
T dr 



727^ ' (3 



= 0. 



(14) 



The adopted solution to the first equation at t <C 1, 
where » 1, is 



, (3 - 2e)i/2 



t 



(15) 



Under the initial condition 3> 1 at r <C 1 the field 
approaches this solution with increasing t. At r ^ 1 
the solution approaches a Bessel function, with am- 
plitude that decays as 



r-3/(2e) 



(16) 



zkx 



0°,|V(2fc). (20) 



Now let us estimate 0^. The wavelength belong- 
ing to comoving wavenumber fc passes through the 
Hubble length at expansion parameter Ofc, where 



(21) 



The second expression follows from the expansion 
law in equation (|l|). The characteristic epoch at 
which ^ passes through unity is ti = 1 (eq. |l5||), 
or t\ = ^^^e^^/^, when i?i = l/(eti). This defines 
the characteristic wavenumber of modes frozen as 4* 
is passing through unity. 



ki/ai = /ie 



-1/2 



(22) 



A mode with wavenumber ~ ki is frozen at the value 
(pl ~ 1/ai (eq |l^), so the rms field fluctuation on 
this scale is frozen at the value 



2.3. The CDM Fluctuation Spectrum 

The CDM field expressed as a quantum field oper- 
ator is 



(x,i) 



E 



(2/cK 



-^</>.(Oe''^- + hc, 



(17) 



where (f)k{t) is the solution to the second of equa- 
tions (p^. The normalization is fixed by the solution 
valid at small r (and for e < 1), where the frequency is 
large compared to the Hubble parameter H (eq. [Q): 



bk — > a{t) ^ exp —i j cudr, 

1 / 



1^^ 



(18) 



cfi ^ ki^'P^.iki) ~ (fci/ai)2 ~ Hi' 



(23) 



Because the field has become effectively massless it is 
frozen at a value on the order of the Hubble param- 
eter at freezing. At shorter wavelengths, the modes 
are frozen later, at smaller H. This is the wanted 
convergence of the field fluctuations on small scales. 

At wavelengths much longer than k^^ the mode 
stops oscillating at ai. In the time interval 

flfc a <C fli the wavenumber in the second part 
of equation (^^ is unimportant and the Higgs field 
varies as ^ oc l/t (eq. ||l^). The result is that the 
mode amplitude varies as a power law of the time. 



(pk oca 

2q = 3-e-[{3-ef- 4e(3 - 2e)^/p] 



1/2 



(24) 
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This is the less rapidly decaying solution to the sec- 
ond of equations ( p^ ) when the term oc may be 
neglected. Thus the amplitude of a long wavelength 
mode is 4>k ~ l/ofc at Hubble crossing (eq. [p^), 
and the amplitude decays by the factor {ak/aiY from 
Hubble crossing to r = 1. After this ^' decreases more 
rapidly than H (eq. |^ ) so the amplitude approaches 
the limiting value 

</)?.~arV«?cxfc-(i-^)/(i-). (25) 

With equations (^ and (p^), the power spectrum is 

= -[(3 - ef - 4e(3 - 2e)^/PY'^/{\ - e), (26) 

on scales much larger than the break at fc fiaic^^^'^ 
(eq. H). 

Figure 1 shows an example of numerical solutions 
for the evolution of the amplitudes of Fourier modes 
of the CDM field. In the computation the amplitude 
is written as 



(j)k{t) — A(t) exp —i / /jjdr, 



(27) 



as in equation (|l8|), and the field equation is writ- 
ten as a differential equation for A{t). The numerical 
computation of A{t) is straightforward. The figure 
shows that \(j>k\'^ decays as l/a{t)'^ until the wave- 
length reaches the Hubble length (eq. [|l8|), then de- 
cays as ^ a"^"? until ^E" starts to decrease rapidly at 
r ~ 1 (eqs. Q and (24j), and then approaches a con- 
stant value, as expected. The mode starts oscillating 
again after inflation when the Hubble parameter has 
decreased to H ^ m, where m is the final CDM mass. 
This is discussed in the next section. 
Figure 2 shows the power law index 



me = dlogP^/dlogfc, 



(28) 



for the power spectrum (eq. [pO[) from the numeri- 
cal solution in Figure 1. The index approaches equa- 
tion ( p6| ) at fc ^ fci , and falls below the scale-invariant 
value me — —3 at fc ^ fci, again as expected. With 
the numerical values in the next section, wavelengths 
of interest for structure formation are some six orders 
of magnitude longer than at the break in the spec- 
trum, and is quite close to a power law with the 
index m^ in equation (p6|). 

The primeval CDM mass distribution is propor- 
tional to the square of 0(x) (eq. ^M). Since 4>{x) is a 



Gaussian process with zero mean value, we have 

(0(xi)20(x2)2) = (02) + 2(0(xi)0(x2))2, (29) 
so the dimensionless mass autocorrelation function is 

(p(xi)p(x2)) ^ (0(X1)0(X2))2 

i;p{xi2) = 1 = 2 — . (30) 



That is, the mass autocorrelation function is propor- 
tional to the square of the field autocorrelation func- 
tion. If the power spectrum of the field fluctuations 
is oc fc"* (eq. H]), with -3/2 > m^ > -3, then 
the Fourier transform is ^0 oc 7'-(3+™«)^ the mass au- 
tocorrelation function is i^p oc oc 7.-(6+2m^)^ ^^^^ 
the Fourier transform is the mass fluctuation power 
spectrum 



Pp oc 



= 3 + 2m^ 



(31) 



Equation (^) indicates the small-scale primeval mass 
fluctuations are mildly nonlinear, ^p(O) = 2 or 5p/ p = 
21/2. 

The fit to the second moments of the large-scale 
mass and CBR distributions requires nip ~ — 1.8 (Pa- 
per II). The integral over this power law diverges at 
large k, so there has to be a break from the power law 
shortward of which Pp decreases more rapidly than 
k^^ . This break is produced by the mass /i of the 
Higgs field, which defines the coherence length of the 
frozen field fluctuations (eq. |0|), 



-iei/2ai-i. 



(32) 



After inflation and at H ^ m the mass density in 
the CDM is determined by the rms value of the field 
frozen at a ai (eq |E3|), 



m2(02)/2 - m20i2 



m2iJi^ 



(33) 



Numerical Values 



This section gives numerical values for the param- 
eters in the inflation model based on a fit to the ob- 
servational constraints. The computations are carried 
to fractional accuracy of first order in e. 

For definiteness I adopt the cosmological model pa- 
rameters 

n = 0.2, A = 0.8, 

To = 2.73 K, Ho = 70 km s^^ Mpc"\ (34) 

and the inflation model parameters 

e = 0.1, 7//3 = 3.3, m^ = -2.4, nip = -1.8. (35) 
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The three famihes of neutrinos are massless. The 
model universe is cosmologicahy flat; A is the contri- 
bution of a term in the stress-energy tensor that acts 
hke Einstein's cosmological constant. The density pa- 
rameter in baryons is JIb, and the density parameter 
in the CDM is O — JIb • The density parameter is re- 
duced from values considered in Peebles (1997) to get 
a better cluster mass function, following the advice 
of Cen (1997). This is discussed further in Paper II. 
The first two parameters in equation (^ produce the 
tilted power law index for the field power spec- 
trum, which translates to the mass fluctuation spec- 
trum power law index rup = — 1.8 (eqs. |^ and [^). 
This gives a reasonable fit to the spectrum of angular 
fluctuations of the thermal cosmic background radi- 
ation and the large-scale fluctuations of the present 
galaxy distribution, which are assumed to be traced 
by the galaxies (Peebles 1997; paper II). The same 
value of 7710 follows from a smaller value of e and a 
larger value for 7//?, a case that is not considered 
here. 

3.1. Parameters for Inflation 

Here I review common measures of the expansion 
factor during and after inflation. If the temperature 
at the end of inflation is T^, and the entropy is shared 
among about as many particles and spin states as 
now, the expansion factor from the end of inflation to 
the present is 



ar 



Tr 
To 



IO^^Tm, 



(36) 



where To = 2.73 K and Tr is expressed in units of 
10" GeV. The mass density - aT^ fixes the Hub- 
ble parameter H^; the ratio to the present value in 
equation ( |35| ) is 



Hr/Ho - lO^^Tl 



(37) 



Fluctuations in (p on the scale of the present Hubble 



length are frozen at epoch a 



when the 



Hubble length is ^. The expansion factor from 
then to the present is 

Oo/flp = Hp/Ho. (38) 

If the power law expansion a oc f^/'^ persists to the 
end of inflation then the expansion factor from ap to 
Qr is 



Qo Ho 



l/(l-e) 



10^^Ti4^\ (39) 



for e = 0.1. This is in the familiar range of values 
(Kolb & Turner 1990). 

3.2. Adopted Inequalities 

The solutions for ip and the CDM fleld (p in §2 are 
approximations based on inequalities to be discussed 
here. The condition 



(40) 



at T < 1 assures that the energy density in (p is small 
compared to that in ?/', and that the coupling to (p has 
little effect on the evolution of V', as assumed in §2. I 
adopt e — 0.1, not a very small number, so the wanted 
power law index requires 7 ~ /?. Since cp is frozen 



at a value comparable to H — 1/ei, equations (15) 
and ( ^ ) require 

e > 7 ~ (41) 

The computation assumes the energy density is 
dominated by that of the field x that is driving in- 
flation (eq. j|]). At T ^ 1 this condition is 



9 ^ 3m2j 



(mpitf > 67r//3. (42) 



The above numbers imply that the time tp of freezing 
of the fluctuations we observe at the present Hubble 
length satisfies 



nipitp 



109 



-24e 



rp2+e 



108Ti4-2-1 > (67r/3-i)i/2, 



(43) 



for e = 0.1. Under the condition to be discussed 
next, that inflation is cool enough that curvature fluc- 
tuations are subdominant, this condition applied at 
the time the density fluctuations we can observe are 
produced is easily reconciled with the condition that 
/3«e (eq. |4l)). 

3.3. Primeval Adiabatic Density Fluctuations 

The isocurvature model assumes primeval curva- 
ture (adiabatic) fluctuations are subdominant on scales 
of interest for structure formation. I consider first the 
effect of the density fluctuations frozen into the CDM 
field 0, and then those of the field x that drives infla- 
tion. 

In time-orthogonal coordinates the components of 
the metric tensor are 500 = 1; ffoa — 0, and 



?Q/3 = -a(i)^ {5af3 - hap) 



(44) 
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In linear perturbation theory the space-time part of 
the Einstein field equation is (Peebles 1980) 

- = SttGToq (p<f>^a/mpi'^. (45) 

To see the meaning of the combination of terms in 
the perturbation hafj to spacetime curvature consider 
the Fourier transform of the field equation. With co- 
ordinates oriented so the wave vector k of a Fourier 
component is parallel to the x'^ axis, equation ( |45| ) is 
ikhii = 87rGro3(k). The component /i33 along the di- 
rection of k does not appear in this equation and the 
value is of no interest because it is freely adjusted by 
a coordinate relabeling along the a:'^ axis. The trans- 
verse part hii — for the inhomogeneous solution 
gives a measure of the curvature fluctuation produced 
by the freezing of the field component, 

hii = -8niG J Toa/kdt - S{t)[a{t)/ktf. (46) 

As discussed in Peebles (1980, §86), and indicated 
in the last expression, hn determines the amplitude 
of the mass density contrast 6 = 6p/p in the most 
rapidly growing adiabatic mode at the end of infla- 
tion. When the proper wavelength again approaches 
the Hubble length, at a/fc ~ t, the observed adiabatic 
density contrast on the scale of the Hubble length is 
{dp/p)h ^ hii. I keep the notation hu for the result 
of summing over Fourier components to get the cur- 
vature fluctuation frozen on the comoving scale k~^, 
which is equal to the adiabatic CBR temperature fluc- 
tuation ST/T hii appearing at the Hubble length. 

As fluctuations in the CDM field (j) are freezing 
the factors (j) ^-nd (t>,a/ci both are on the order of 
Hcj} ~ H^, so equation ( p5[ ) indicates the tempera- 
ture fluctuation appearing now at the Hubble length 
is 

hn {Hp/mpif ^ 10-2"+54^Ti4^ « 10"^ (47) 

The last expression is the condition that the adia- 
batic perturbation resulting from the freezing of (j) is 
small compared to the isocurvature perturbation that 
will be tuned to match the observations. The same 
condition eliminates the tensor quantum fluctuations 
frozen into the spacetime curvature. 

The adiabatic perturbation from the fluctuations 
frozen into the fleld that is driving the power law in- 
flation produces CBR temperature anisotropy 

ST/T - {Tr/nipif < 10-^ (48) 



At e = 0.1 and T14 1 the conditions in equa- 
tions (^^ and ( ^ ) are well satisfied; curvature fiuc- 
tuations are negligibly small. 

3.4. The Coherence Length of the CDM Field 

The parameters for the inflation model are con- 
strained by the coherence length of 4>{x) at the break 
in its power law spectrum (eq. [^): the coher- 
ence length flxes the mean square value of the fleld, 
which in turn fixes the mean mass density in CDM, 

P = TO^(0^). 

The required value of the coherence length follows 
from the normalization to the mass fiuctuation spec- 
trum. The spectrum extrapolated to the present in 
linear perturbation theory is taken to be 

Ppik) = J dV^pe"'-- 

= 6300h~^{0.1h/ky-^ Mpc^ (49) 

This is consistent with the assumption that galaxies 
trace mass on large scales (Paper II). The unit for 
the wavenumber is radians per megaparsec, and the 
power law index is the adopted value in equation ( ^5| ) . 
The Fourier transform of Pp is the mass autocorrela- 
tion function in linear perturbation theory, 

= 7.2(/irMpc)"''. (50) 

The primeval mass distribution is homogeneous, 
with fiuctuations in the local CDM mass density bal- 
anced by opposite fiuctuations in the baryons and 
radiation. As these fluctuations in composition ap- 
pear at the Hubble length the radiation distribution 
is smoothed, at short wavelengths by the propagation 
of acoustic waves that are dissipated by photon diffu- 
sion in the baryon-radiation fluid before decoupling, 
and at large wavelengths by free streaming of the ra- 
diation after decoupling. This leaves the CDM as 
an isothermal distribution within the Hubble length. 
Wavelengths in the CDM distribution that reach the 
Hubble length prior to the epoch Zeq of equality of 
mass densities of matter and radiation are left with 
close to the primeval amplitude at z = z^q. At longer 
wavelengths the CDM distribution at Zcq is smoother 
than the initial value because the net mass density re- 
mains close to homogeneous when pressure gradient 
forces may be neglected. 

At wavelengths small compared to the Hubble 
length at z = Zcq the amplitude of a CDM fluctu- 
ation varies with time as the isothermal perturbation 
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equation 

D + -D^A7:GDpcDM = U-Y^DHo', (51) 
a 2 \ a y 

where the expansion rate equation is 
(a/a)2 = Ho'' 



n {ao/af' + il,. (flo/a)^ + A 



(52) 



The density parameter in radiation is 
Vlr = 4.2 X IQ-^h-'^. 



(53) 



The density parameter in baryons plus CDM is fi, and 
equation (|5^) assumes the CDM mass density is much 
larger than that of the baryons. A more complete 
analysis would describe the collapse of the baryons 
into the CDM mass concentrations after decoupling, 
but the resulting correction to equation ( 51 ) is beyond 
the level of accuracy of this discussion. Equations ( |5l| ) 
and ( p2| ) assume a cosmologically flat model, with A = 

Results of numerical integration of these equations 
with h = 0.7 are shown in Figure 3. The net growth 
factor from high redshift to the present is well approx- 
imated by the fitting function 

e{n) = D{z oo)/D{z = 0) = 5.7 x 10"^r2"^-^^ 

(54) 

In linear perturbation theory the primeval CDM mass 
autocorrelation function is e(f2)^ times the present 
value in equation ( |50| ) at fixed comoving length. The 
field autocorrelation function is related to the mass 
function by equation (^0|). The result is the primeval 
field autocorrelation function. 



2) 



1 X 10" 



i^) (/iXMpc)°-6r!l-22- 

At f2 = 0.2 this power law extrapolates to unity at 



(55) 



hxi — 7 pc. 



(56) 



In the approximation that the shape of the power 
spectrum changes fairly rapidly from Pp{k) oc 
on large scales to Pp oc with s < —3 on smaller 
scales, equation ( |56|) is the required comoving coher- 
ence length of the CDM field 0(x). 

3.5. The CDM Mass and Mass Density 

Now let us consider the condition that the energy 
density in left from inflation matches the present 
mean mass density in the cosmological model. 



The present cosmological mass density is 

Po = 1.88 X IQ-'^^nh' g cm-3 (57) 
= 8 X IQ-^'^nh' GeV^ (58) 
= m^^o^, (59) 

where m is the field mass and 4>o is the present rms 
value of the field. This gives 



bo = 10"23;j^l/2^g-l QgY^ 



(60) 



where 

m = niQ Gev. (61) 

The field fiuctuations are nonrelativistic, because the 
de Broglie wavelength at the coherence length a{t)xi 
is much longer than the Compton wavelength all the 
way back to inflation, so at < to the rms field 
value varies as </> oc a(t)^'^^'. When the universe is 
radiation-dominated the Hubble parameter is 

iJ - 2 X 10-'°{ao/af s"^ 

= 1 X I0-'^^{ao/af GeV, (62) 

if the number of particle families and spins at very 
high redshift is about the same as at z ~ 10^*^. The 
field starts oscillating after inflation at ~ to, when 
the expansion factor is 

Oo/a™ ~ 1 X lO^^TOgl/^. (63) 

The field thus is frozen at the rms value 

Hi - (1)1 (t>o{ao/a,nf''' 

~ 7 X lOVg-i/^/il^i/^ (64) 
~ 2 X lOVg-i/'' GeV, 

for the cosmological parameters in equation (^4|). 

Equation (^) fixes the value of the Hubble param- 
eter when the field coherence length is equal to the 
Hubble length (eq. [^). The model assumes power 
law expansion, H oc a~'^, continues to the end of infla- 
tion at expansion parameter a^, so the present value 
of the coherence length is 



ai 



(65) 



With equations (|3^) and (^ and e = 0.1 this gives 
the expansion factor from freezing of the field at its 
coherence length to the end of inflation. 



ar/ai ~ 10i6[ri(pc)Ti4]i-^ 



(66) 
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This is much smaller than the expansion factor at 
freezing of the field fluctuations on the scale of the 
present Hubble length, of course (eq [p9[). 

The rms value of the frozen field is (pi ^ Hi ^ 
Hr{ar / aiY . The temperature at the end of inflation 
follows from equations (|3^), (^), and (|6^): 

Ti4-^0.2m9'°iVi(pc)~"05. (67) 

This indicates that for a considerable range of values 
of the CDM particle mass, m, inflation would have to 
end at a relatively cool temperature, about 10^'^ GeV. 
Finally, the mass /i in the Lagrangian (|^) is given by 
equations ( |3^ ) and ([6^), 

^ - e^/^Hi - lOVg-i/^ GeV (68) 

4. Discussion 

In the adiabatic CDM family of models (that I have 
termed ACDM) the cold dark matter particles can 
have been produced out of entropy at any time af- 
ter inflation ended, a broadly general situation. The 
isocurvature CDM (ICDM) model is more specific — 
the CDM is a remnant of a field present during infla- 
tion — and hence arguably less likely in the absence 
of evidence for or against the nature of the CDM. This 
consideration may become more compelling when we 
have better experimental constraints on the dark mat- 
ter. 

Assessments of the relative degree of "fine tuning" 
of the adiabatic and isocurvature inflation models also 
seem to be of doubtful significance at this stage of 
our understanding, so I note only four aspects of the 
ICDM model. First, there is precedent in the litera- 
ture for the CDM potential 0) in equation (|l0|). 
The same functional form, applied in a different con- 
text, appears in the Hybrid Inflation model (Copeland 
et al. 1994; Linde 1994), and similar forms are not 
uncommon (eg. Kofman & Linde 1987; Hodges et al. 
1990; Randall, Soljacic & Guth 1996). Second, the 
conditions on an acceptable set of dimensionless pa- 
rameters for the ICDM model are 7 ~ /? ^ e ~ 0.1. 
These do not seem unduly severe or artificial. Third, 
the model requires that the initial value of the field 
ip be large enough that approaches the solution in 
equation ( p^ ) well before field fluctuations are frozen 
on scales we can observe, and that the initial value of 
(p be small enough that ip can drive it to zero mean 
value before the field is squeezed on scales of interest. 
I do not know how to judge whether these conditions 



are likely outcomes of the physical situation prior to 
inflation. Finally, an acceptable set of values of the 
characteristic masses that appear in the model is 

M - 10^^ GeV, Tr - 10" GeV, 

/I - 10^ GeV, m - 1 GeV. (69) 

The broad range of values is impressive, but so is the 
range of measured characteristic masses in particle 
physics and those that appear in physics-motivated 
models for inflation (Randall 1997). 

The greatly tightened constraints from observa- 
tional programs in progress are going to make it much 
harder to invent acceptable structure formation mod- 
els such as ACDM and ICDM from largely theoreti- 
cal and aesthetic considerations. Some already known 
member of the CDM family — in which I would in- 
clude ACDM and ICDM — may survive the precision 
tests in progress and by its success compel acceptance 
as a truly valid approximation to reality. If this does 
not happen, perhaps we will be lucky enough to see 
in the results from observations and particle physics 
guidance to the formulation of more promising mod- 
els. 

I conclude that the ICDM model has as valid a 
pedigree from accepted ideas about the very early 
universe as might be expected for a model that is mo- 
tivated by astronomical phenomenology rather than 
physics. The observational situation is discussed in 
Paper II. 

This work was supported in part by the NSF. 
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Fig. 1. — Evolution of amplitudes of Fourier modes 
of the CDM field for e = 0.1 and 7 = 3.3/3. The 
dimensionless wavcnumbcrs are k = 10^**, 10^^, 1, 
and 100 from top to bottom. 
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Fig. 2. — Logarithmic derivative of the CDM field 
power spectrum (eq. p^) for a numerical solution 
with the parameters in Figure |l|. This shows the re- 
gion near the break in the spectrum. 
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Fig. 3. — Evolution of the CDM density contrast 
on scales small compared to the Hubble length at 
redshift Zeq, in linear perturbation theory. The 
primeval isocurvature density fluctuation becomes an 
isothermal CDM fluctuation when its length scale be- 
comes comparable to the Hubble length and radiation 
pressure smooths the distribution of radiation and 
baryons. The isothermal density fluctuation starts 
to grow as the redshift approaches Zeq. 
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